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We analyze the isospin-breaking corrections to quark condensates within one-loop SU{2) and 
5*17(3) Chiral Perturbation Theory including m„ 7^ md as well as electromagnetic (EM) contribu- 
tions. The explicit expressions are given and several phenomenological aspects are studied. We 
analyze the sensitivity of recent condensate determinations to the EM low-energy constants (LEG). 
If the explicit chiral symmetry breaking induced by EM terms generates a ferromagnetic-like re- 
sponse of the vacuum, as in the case of quark masses, the increasing of the order parameter implies 
constraints for the EM LEG, which we check with different estimates in the literature. In addition, 
we extend the sum rule relating quark condensate ratios in SU (3) to include EM corrections, which 
are of the same order as the m„ 7^ md ones, and we use that sum rule to estimate the vacuum 
asymmetry within GhPT. We also discuss the matching conditions between the SU{2) and SU{3) 
LEG involved in the condensates, when both isospin-breaking sources are taken into account. 
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The low-energy sector of QCD has been successfully described over recent years within the chiral lagrangian 
framework. Chiral Perturbation Theory (ChPT) is based on the spontaneous breaking of chiral symmetry 
SUL{Nf) X SUR{Nf) — SUviNf) with Nf = 2,3 light flavours and provides a consistent and systematic model- 
independent scheme to calculate low-energy observables [IH3]- The effective ChPT lagrangian is constructed as the 
more general expansion C = Cp2 + Cp4 + . . . compatible with the QCD underlying symmetries, where p denotes 
derivatives or meson mass and external momentum below the chiral scale A,^, ^ 1 GeV. 

The SUv{Nf) group of vector transformations corresponds to the isospin symmetry for Nf = 2. In the Nf = 3 
case, the vector group symmetry is broken by the strange-light quark mass difference — m„^d, although nis can 
still be considered as a perturbation compared to A^, leading to SU{3) ChPT [3j. In the Nf — 2 case, the isospin 
symmetric limit is a very good approximation in Nature. However, there are several known examples where isospin 
breaking is phenomenologically relevant at low energies, such as sum rules for quark condensates |3j, meson masses 
and corrections to Dashen's theorem [Ij , pion-pion [3 IS] and pion-kaon [71 [S] scattering in connection with mesonic 
atoms [H [To] , CP violation ]TT] , ap — /o mixing [T2] , kaon decays [l3|, [14] and other hadronic observables (see [H] for 
a recent review). 

The two possible sources of isospin breaking are the nid — to„ light quark mass difference and electromagnetic (EM) 
interactions. Both can be accommodated within the ChPT framework. The former is accounted for by modifying the 
quark mass matrix and generates a tt^t] mixing term in the SU{3) lagrangian pj. The expected corrections from this 
source are of order (m^j — m„)/ms. On the other hand, EM interactions, which in particular induce mass differences 
between charged and neutral light mesons, can be included in ChPT via the external source method and give rise to 
new terms in the effective lagrangian [4H61 IT6HI8] of order £^2 , Ce^pi and so on, with e the electric charge. These terms 
fit into the ChPT power counting scheme by considering formally = 0{p'^ /F'^), with F the pion decay constant in 
the chiral limit. 

The purpose of this paper is to study the isospin-breaking corrections to quark condensates, whose main importance 
is their relation to the symmetry properties of the QCD vacuum. The singlet contributions {uu -f dd) for SU (2) and 
{uu + dd + ss) for SU(3) are order parameters for chiral symmetry, while the isovector one {uu — dd) behaves as 
an order parameter for isospin breaking, which is not spontaneously broken [TS]. We will calculate the condensates 
within one-loop ChPT, which ensures the model independency of our results, and will address several phenomenological 
consequences. The two sources of isospin-breaking will be treated consistently on the same footing, which will allow us 
to test the sensibility of previous phenomenological analysis to the EM low-energy constants (LEC). Moreover, the EM 
corrections induce an explicit breaking of chiral symmetry which will lead to lower bounds for certain combinations 
of the LEC involved, provided the vacuum response is ferromagnetic, as in the case of quark masses. In addition, in 
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SU{3) one can derive a sum rule relating the different condensate ratios for niu rud [5] which, as we will show here, 
receives an EM correction not considered before and of the same order as that proportional to m„ — m^. The latter is 
useful to estimate the vacuum asymmetry {dd) / (uu) reliably within ChPT. An additional aspect that we will discuss 
is the matching of the SU{2) and SU{3) LEG combinations appearing in the condensates when both isospin-breaking 
sources are present, comparing with previous results in the literature. The analysis carried out in the present work 
will serve also to establish a firm phenomenological basis for its extension to finite temperature, in order to study 
different aspects related to chiral symmetry restoration [20 . 

With the above motivations in mind, the paper is organized as follows: in section [IT] we briefly review the effective 
lagrangian formalism needed for our present work, paying special attention to several theoretical issues and to the 
numerical values of the parameters and LEG needed here. Quark condensates for SU{2) are calculated and analyzed 
in section [Tin where we discuss the general aspects of the bounds for the EM LEG based on chiral symmetry breaking. 
In that section we also comment on the analogy with lattice analysis. The SU (3) case is separately studied in section 



IV In that section, we perform first a numerical analysis of the isospin-breaking corrections, paying special attention 
to the effect of the EM LEG in connection with previous results in the literature. In addition, we obtain the EM 
corrections to the sum rule for condensate ratios, which we use to estimate the vacuum asymmetry within GhPT. 
We also provide the LEG bounds for this case, checking them with previous LEG estimates and, finally, we discuss 
the matching conditions for the LEG involved. In Appendix |A] we collect the lagrangians of fourth order and the 
renormalization of the LEG used in the main text. 

II. FORMALISM: EFFECTIVE LAGRANGIANS FOR ISOSPIN BREAKING 

The effective chiral lagrangian up to fourth order is given schematically by 

^eff = ^p-^+e^ + ^p^+e^p^+e^- (1) 

The second order lagrangian is the familiar non-linear sigma model, including now the gauge coupling of mesons to 
the EM field through the covariant derivative, plus an extra term: 

£p2+,2 = — tr [D^U^D^'U + 2B^M {U + U^)] 
+ Ctr [QUQU^ . 

Here, F is the pion decay constant in the chiral limit and U{x) G SU{Nf) is the Goldstone Boson (GB) field in the 
exponential representation U = exp[i$/_F] with : 



SU{2) : $ 



\V2tt- -TfO 



j^rj V2t:+ V2K^ 



SU{3) : $ = V2Tr- -n" + -j=r^ V2X" , (3) 



V2K- V2K° 3|'7 



V3' 

7S' 



with 77 the octet member with = S ~ 0. The covariant derivative is Z?^ = 9^ + iAf^[Q, ■] with A the EM field. 
M and Q are the quark mass and charge matrices, i.e., in SU{3) Ai = diag(mti, m^, m^) and Q — diag(ei,, Cd, Cg) 
with Cu = 2e/3, — Cs ~ — e/3 for physical quarks. The additional term in ([2]), the one proportional to C, can be 
understood as follows: the QGD lagrangian for m„ = coupled to the EM field is not invariant under an isospin 
transformation q — !■ gq with g e SU{Nf) and q the quark field. However, it would be isospin invariant if the quark 
matrix Q is treated as an external field transforming as Q — > g^Qg- Therefore, the low-energy effective lagrangian 
has to include all possible terms compatible with this new symmetry, in addition to the standard QGD symmetries. 
The lowest order 0{e^) is the C-term in since U transforms as [/ — >■ g^Ug. Actually, one allows for independent 
"spurion" fields Ql{x) and Qr{x) transforming under SUL{Nf) x 5'C/_r(7V/) so that one can build up the new possible 
terms to any order in the chiral lagrangian expansion, taking in the end Ql = Qb = Q 

In the previous expressions, F, Bomu^d.s,C are the low-energy parameters to this order. Working out the kinetic 
terms, they can be directly related to the leading-order tree level values for the decay constants and masses of the 



3 



GB. In SU{2)^ the tree level masses to leading order are: 

Ml^ = Ml- ^2mBo + 2C^, 

Mlo = 2mBo, (4) 

with m = (to„ + md)l2 the average light quark mass. Note that both terms contributing to the charged pion mass 
are of the same order in the chiral power counting, although numerically (M^± — Af^o) /A^^o — 0.1, which we will 
use in practice as a further perturbative parameter to simplify some of the results. 

In the S'[/(3) case, the mass term in ([2| induces a mixing contribution between the 7r° and the i] meson fields 
given by Cmix — [B^ / y/2>){'md — 'mu)'n'^rj. Therefore, the kinetic term has to be brought to the canonical form before 
identifying the GB masses, which is performed by the field rotation [3]: 

7r° = 7f° cos £ — ?7 sine, 
rj = 7f° sine + 77 cose, (5) 

where the mixing angle is given by: 



tan 2e = 



\/3 mrf 



2 TTis — m 

Once the above n^rj rotation is carried out, the SU{i) tree level meson masses to leading order read: 



(6) 



= Ml- 
= 2Bo 



2mBo + 2C 



2 , sin^ e 

m m , —ml 

3^ '^cos2e 



= Ml_ = (m, + m„)So + 2C 



rj 



(nis + md)Bo 
2Bo 



1/^ „ N 2, 
-(m + 2m, H — (m. 



p2 ' 



m) 



sin^ e 
cos 2e 



(7) 



The above five equations are the extension of the Gell-Mann-Oakes-Renner (GOR) relations [21] to the isospin 
asymmetric case and allow to relate the four constants Bomu,d,s and C (e is given in terms of quark masses in ^) 
with the five meson masses or combinations of them. The additional equation provides the following relation between 
the tree level LO masses: 



(a4± - Ml^Y - 3 (Af2 _ M|„) (M|,o - Mlo) = 0, (8) 

The above equation is compatible with the one obtained in neglecting 0{mu — rrid)'^ terms. Actually, note that 
although all terms in ([t]) are formally of the same chiral order, numerically (see below) we expect e ~ {Vi/4){md — 
m,()/ms <C 1 and hence the mixing-angle corrections to the squared masses to be 0{Mle) and ©(Af^e^) for the 
neutral pion and eta respectively. On the other hand, in the isospin-symmetric limit (m„ = m^ and e = 0) (|8]) is 
nothing but the Gell-Mann-Okubo formula 4Af|. - 3M^ - Af^ = Q [23]. Neglecting only the md — m„ mass difference 
in ^ leads to Dashen's theorem A/|.± - M]^o = Af^i - Af^o [24] and then eq. (|8| reduces to 4Ar|,o - 3Af2 _ = g, 
i.e.Ixhe Gell-Mann-Okubo formula for neutral states. However, the violation oiDashen's theorem at tree level due to 
those quark mass differences is significant numerically for kaons. In our present treatment we consider those differences 
on the same footing as the EM corrections to the masses. For pions, the main effect in the n'^ — 7r+ mass difference 
comes from the EM contribution |25j . 

All the previous expressions hold for tree level LO masses M^ with a = tt* , tt", , ry, in terms of which we will write 
all of our results. They coincide with the physical masses to leading order in ChPT, i.e., M^ ^f^^^ = M^{1 + 0(Af^)). 
Calculating the GhPT corrections to a given order allows then to determine the numerical values of the tree level 
masses, knowing their physical values and to that order of approximation. The same holds for F, which coincides with 
the meson decay constants in the chiral limit F^ ,^ — F^{1 + ©(Af^)). Next to leading order O(Af^) corrections to 
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meson masses and decay constants were given in [51 [3] for = 0. EM corrections to the masses can be found in [3] 
for SU{3) and in [6l[T8] for SU{2) including both ^ and m„ ^ nid isospin-breaking terms. 

The fourth-order lagrangian in ([T]) consists of aU possible terms compatible with the QCD symmetries to that order, 
including the EM ones. The Cpi lagrangian is given in for the SU{2) case, /ii.2,3 (contact terms) and /1...7 denoting 
the dimensionless LEG multiplying each independent term, and in 3 for SU{3) the LEG named H12 and £i...io. 
The EM £e2p2 and for SU{2) are given in OlSj, fci,...i3 denoting the corresponding EM LEG, and in \i. for SU{3) 
with the Kx n EM LEG. For completeness, we give in Appendix \K\ the relevant terms needed in this work. The LEG 
are renormalized in such a way that they absorb all the one-loop ultraviolet divergences coming from Cp2 and C^i , 
according to the GhPT counting, and depend on the MS low-energy renormalization scale in such a way that the 
physical quantities are finite and scale- independent. The renormalization conditions for all the LEG can be found in 
[H El ini HZ] and we collect in Appendix [A] only those needed in the present work. 

As customary, we denote the scale-dependent and renormalized LEG by an "r" superscript. The renormalized 
LEG are independent of the quark masses by definition, although their finite parts are unknown, i.e., they are not 
provided within the low-energy theory. The numerical values of the LEG at a given scale can be estimated by fitting 
meson experimental data, theoretically by matching the underlying theory under some approximations, or from the 
lattice. These procedures allow to obtain estimates for the "real-world" LEG at the expense of introducing residual 
dependencies of those LEG on the parameters of the approximation procedure, which typically involves a truncation 
of some kind. Examples of these are the nis dependence on the SU{2) LEG when matching the SU{3) ones, the 
correlations between LEG, masses and decay constants through the fitting procedure, the QGD renormalization scale 
and gauge dependence of some of the EM LEG or the dependence with lattice artifacts such as finite size or spurious 
meson masses. We will give more details below, specially regarding the EM LEG which will play an important 
role in our present work. An exception to the LEG estimates are the contact LEG hi and Hi, which are needed for 
renormalization but cannot be directly measured. The physical quantities depending on them are therefore ambiguous, 
which comes from the definition of the condensates in QGD perturbation theory, requiring subtractions to converge 
[2]. It is therefore phenomenologically convenient to define suitable combinations which are independent of the hi, Hi. 
We will bear this in mind throughout this work, providing such combinations when isospin-breaking is included. 

We will analyze in one-loop GhPT (next to leading order) the quark condensates, which for a given flavour Qi can 
be written at that order as: 



ff 



drrii 



(9) 



The above equation is nothing but the functional derivative with respect to the i-th component of the scalar current, 
particularized to the values of the physical quark masses, according to the external source method [21 [3] ■ Therefore, 
we will be interested only in the terms of the fourth-order lagrangian containing at least one power of the quark 
masses. These are the operators given in eqns.( [Al] ) and (A2| for SU{2) and SU{?>) respectively. Thus, the LEG that 
enter in our calculation are Z3, ft-i, /i3, fcs, /cg? ^7 in SU{2), and L^, Lg, H2, K^, K^, Kq, Kiq in SU{3). Besides, up to 
NLO, only tree level diagrams from the fourth-order lagrangian can contribute to the condensates, so that in practice 
it is enough to set f/ = 1 in ( Al )-( A2 ) for getting those tree- level contributions from 



A. Masses and low-energy constants 



For most of the numerical values of the different low-energy constants and parameters in the SU{5) case, we will 
follow [55], where fits to K14 experimental data are performed in terms of 0{p^) GhPT expressions, including the 
isospin mass difference m„/md ^ 1 and EM corrections to the meson masses, extending a previous work [27j where 
isospin breaking was not considered. Those fits have been improved in a recent work ^28) . which takes into account 
new phenomenological and lattice results. We will however stick to the values of [26], since our main interest is to 
compare the isospin-breaking condensates with the two sources included and to estimate the effect of the EM LEG. 
In the new fits (28[ isospin breaking is included only to correct for the charged kaon mass and the condensate values 
are not provided. For a review of different estimates of the quark masses and condensates see also [21] and [SU]. In 
addition, in [31' a recent update of lattice results for low-energy parameters can be found, including LEG and the 
quark condensate. We will use the central values of the main fit in ^26J. The value of nig/rh — 24 j29j 30] is used as 
an input in [55], as well as Lg = 0, as follows e.g. from OZI rule or large- A^c arguments [3]. The more recent fits [2S] 
consider an updated value of rris/rh = 27.8 in accordance with recent determinations [311 and a nonzero value of Lg is 
obtained as an output. The suppression of Lg has been questioned in connection with a reduction of the light quark 
condensate when the number of flavours is increased |32j . within the framework of generalized GhPT. In that context, 
the chiral power counting is modifled due to the smallness of the condensate. Here, we will adhere to the standard 
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ChPT picture, where the condensate and the GOR-hke relations are dominated by the leading order [53], sustained 
by the recent lattice LEG estimates [3l!- The values of = 87.1 MeV, 2Bom = 0.0136 GeV^, m„/md = 0.46 and 
Lsin = 770 MeV) = 0.62 x lO'^ are outputs from the main fit in [15]. With those values we get from ^ e ^ 0.014 
and from Q the tree level masses of tt", fC", rj. 

To calculate the tree level charged meson masses, we need also the value of the C constant, which can be inferred 
also from the results in [26 since the EM correction is numerically very small in the charged kaon mass with respect 
to the pure QGD contribution. This allows to extract the tree level charged kaon mass directly from the expressions 
for Mj^± /Mj^± QQD in |26j . approximating Mj^± qqd by the full physical mass. From there we extract the value of 
C by subtracting the tree level QGD part in ^ calculated with the above given quark masses. In this way we get 
C — 5.84 X 10^ MeV*, which is very close to the values obtained simply from the charged-neutral pion mass difference 
in ([t]) setting the masses and F to their physical values [5] or from resonance saturation arguments . From that C 
value we obtain the tree-level charged pion mass, using again ([t]). Nevertheless, to the order we are calculating we 
could have used as well the physical meson masses and decay constants instead of the tree level ones, since formally 
the difference is hidden in higher orders. The main reason why we choose the values in |26j is to compare directly with 
their numerical quark condensates and estimate the importance of the corrections (see section IV A for details) . 
The constant H2 will also appear explicitly in quark condensates. Since it cannot be fixed with meson experimental 
data, when needed we will estimate it from scalar resonance saturation arguments as iJJ = 2Lg |16l 126) , although 
we will comment below more about the H2 dependence of the results and provide physical quantities which are 
independent of the contact terms. 

Regarding the EM LEG, the SU{3) have been estimated in the literature under different theoretical schemes. 
Resonance saturation was used in [35], large- TVc and NJL models in [34], complemented with QGD perturbative 
information in [36j and a sum-rule approach combined with low-lying resonance saturation has been followed in 
[37K35]. The works [Ml f55H55] have in common the use of perturbative QGD methods for the short-distance part of 
the LEG and different model approaches for the long-distance part. This procedure implies that the LEG estimated in 
that way depend (roughly logarithmically) in general on the QGD renormalization scale, which we call /xq to distinguish 
it from the low-energy scale /i, as well as on the gauge parameter. A closely related problem is that the separation 
of the strong (e — 0) and EM contributions for a given physical quantity is in principle ambiguous [331 ISZl |321 HD] . 
The origin of this ambiguity [ID] is that QGD scaling quantities such as quark masses contain also EM contributions 
through the Renormalization Group (RG) evolution in the full QGD-I-EM theory. Thus, a particular prescription 
for disentangling those contributions must be provided. In addition, when matching such quantities between the 
low-energy sector and the underlying theory, the choice of a given prescription will necessarily affect the scale and 
gauge dependence of the EM LEG. These theoretical uncertainties, as well numerical ones, make those theoretical 
EM LEG estimates not fully compatible among them. For these reasons, in many works analyzing EM corrections, 
the LEG are simply assumed to lie within "natural" values \Kl\, — 6.3 x 10^^ at the scale fi ^ Mp [HIS]. 

The above theoretical issues will be addressed in more detail in sections |lll| and [IV D| The LEG dependencies on the 
QGD scale and on the gauge parameter do no affect directly our results, being only relevant when comparing them 
with approaches where those LEG are obtained by matching the underlying theory. In that context we will see that 
the LEG combinations that we will deal with are gauge independent and lie within the stability range where the 
dependence on fiQ is smooth and the matching makes sense |341l36j . The theoretical errors quoted e.g. in |34) account 
for the uncertainty related to the fiQ dependence. 

As for the SU{2) no direct estimate is available to our knowledge, although one can relate them to the 
by performing formally a l/rris expansion in a given physical quantity calculated in SU{3) and comparing to the 
corresponding SU{2) expression, similarly as the k O Li conversion given in [3]. This has been done partially for 
some combinations of the LEG, namely, those appearing in the neutral pion mass |41j . in pion scattering [9l|41] and in 
the pionium lifetime [41, 42j. More recently, a full matching of the EM SU{2) and SU{3) LEG at the lagrangian level 
has been performed in ^43] using functional integral methods in the chiral limit. In the present work, we will provide 
a complementary analysis. Namely, in section IV we will obtain the matching relations between the LEG involved in 
the quark condensates, including both to„ — rud and EM contributions. Those relations will be consistent with the 
results in [43 and phenomenologically useful when dealing with approximate LEG determinations where isospin and 
mass corrections may be entangled. 



III. TWO-FLAVOUR QUARK CONDENSATES AND BOUNDS FOR THE EM LEG 



We start by giving the explicit one-loop GhPT expressions for the quark condensates in SU{2) with all the isospin- 
breaking corrections included, which we derive from ([9]): 
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dd) = -2F^Bq 



][f2 



(uu-dd) = AB^imd-m^)h3- -F^Boe^k7 + 0{p^) , 



(10) 
(11) 



where: 



(12) 



and throughout this work we wih follow the same notation as in [21 [3] : 



327r2F2 



log 



Ml 



327r2 



1 



log^ 



(13) 



The /ii arise from the finite part of the one-loop tadpole- like contribution Gi{x — 0), with G the free meson 
propagator [5]. The renormalization conditions for the LEG involved in (lOl-(ll) can be found in [2 and [6 and we 
give them in Appendix [A| With that LEG renormalization, one can check that the condensates in (lOl-(ll) are finite 
and scale-independent, which is a nontrivial consistency check. Recall that and k^ do not need to be renormalized. 
The condensates still depend on the hi and contact LEG, which, as explained above, yield an ambiguity in the 
determination of the condensates. The result ( 10 1 for e 7^ and 7^ nid reduces for e 
in [2]. The condensate difference (11) is given in ,2J for e 



to the expressions given 
0, which we agree with, and in [6 for e 7^ 0, which we 
also agree with, except for the relative sign between the two terms, which should be a minus in their eq.(3.12) |44j . 

Let us now argue on how the EM corrections to the condensates may lead to constraints for the EM LEG. Those 
corrections come directly from the coupling of the EM field to the quarks in the QGD action, which break chiral 
symmetry. Actually, to understand better the origin of the different sources involved in chiral symmetry and isospin 
breaking, it is useful to keep the charges of the u and d flavours arbitrary and to separate the isoscalar and isovector 
contributions of the charge matrix in SU(2): 



Q 



1 ■ 



ed 



-T3, 



(14) 



with Ta — diag(l, —1) corresponding to the third isospin component. The EM part of the QGD lagrangian q'y^^A^Qq 
breaks explicitly chiral symmetry SUl{2) x SUii{2) if e„ 7^ ed, through the isovector part in (14|. The isovector also 



breaks the isospin symmetry SUv{2) (L — R) except for transformations in the third direction, which corresponds to 
electric charge conservation. On the other hand, the mass term qA4q breaks chiral symmetry for any nonzero value of 
the quark masses, preserving isospin symmetry if m„ = rud- Altogether, the conclusion is that the QGD lagrangian 
is chiral invariant only if e„ — Cd and rriu = rud = 0. Thus, chiral symmetry is explicitly broken even if e„ = e^, as 
long as any of the quark masses 7^ 0, or equivalently, in the presence of an external scalar source, as needed to 
derive the condensates. If e„ = Cd and ruu = rud 7^ 0, chiral symmetry is broken but isospin symmetry is conserved. 

Now, let us remind how this charge and mass symmetry breaking pattern translates into the low-energy sector. 
The leading order £p2_^.^2 in ([2]) contains separate combinations of the charge and mass terms, both sharing the 
QGD pattern. Thus, the charge contribution proportional to C in ([2]) can be decomposed according to ([m]), giving a 
constant term proportional to (e„ + e(j)^ independent of masses and fields, plus the term C [(e„ — 6^)^/4] trjTsUTsW] , 
which contributes directly to the pion EM mass difference in Q. Therefore, in the second order lagrangian all the 
EM chiral symmetry breaking terms are proportional to (e„ — Cd)"^- This is no longer true for the fourth order 
lagrangian in ( Al I, for which the symmetries of the theory allow for crossed mass-charge terms, like those proportional 
to k^jk^ and fcy. Those crossed terms break chiral symmetry even for = Cd for any nonzero quark mass, the 
strength of chiral breaking being proportional both to the quark charge e and to the quark mass rfi. Gonsequently, 
they contribute to {qq) = — 2(9£/9(m„ + nid)), the expectation value of the SUv{2) singlet behaving as an order 
parameter for chiral symmetry breaking. Setting U — t in the lagrangian gives a piece proportional to k^ -l-fcg yielding 
(gu -I- ed)2 and (e„ — e^)^ contributions to (qq) and another one proportional to kr contributing both to (qq) and to 
{uu — dd) = —2{dC/d{mu — fnd)), the isotriplet order parameter of isospin breaking: 



= 2F'^Bo {(fcs -t- ke) [{cu + Cdf + (e„ - Cd^] + nid) 

+ 2kj [(e„ -I- edYirUu + rud) + (e„ - ed)(e„ -|- ed){mu - m^)] } -|- 
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where the dots indicate terms not contributing to the condensates at this order. 

Thus, we see how the two sources of isospin breaking show up in the order parameter (11), which does not receive 
pion loop contributions in SU{2). The latter is the explicit confirmation that isospin symmetry is not spontaneously 
broken in QCD [19J, since all the contributions to this order parameter vanish for m.^ = irid and e = 0. 

On the other hand, the quark condensate (qq) behaves as an order parameter for chiral symmetry and therefore 
measures the different sources of symmetry breaking; spontaneous and explicit. Thus, it is naturally expected that 
its absolute value increases when a new symmetry-breaking source, as the EM contribution, is switched on. This 
increasing behaviour is what we denote as ferromagnetic, in analogy with the behaviour of the magnetization in a 
ferromagnetic material under an external magnetic field. 

There is no a priori formal argument to ensure the ferromagnetic-like nature of the QCD low-energy vacuum. We 
can nevertheless learn from the response of the system to the light quark mass m, which is the actual counterpart of 
the magnetic field in a ferromagnet, since it breaks the chiral symmetry explicitly by coupling to the order parameter. 
For the mass, this ferromagnetic behaviour is actually followed by the e = condensates in the standard ChPT 
framework, both to 0{p^) and to ©(p^) with the LEG in [27], although assuming Lg suppression and the dependence 
on contact terms still introducing a source of ambiguity. With the recent fit giving nonzero Lg and a new value 
also for Lg [5H] we get also that {uu + dd) increases with light quark masses, from our 0{p^) expressions. The same 
ferromagnetic effect implies the increasing of the critical temperature of chiral restoration when increasing the pion 
mass, confirmed by ChPT calculations [45] and lattice simulations [46|. Finally, lattice results for the condensate 
[5T] reveal a systematic increase of its absolute value with respect to direct estimates reflecting again the same 
behaviour, since the pion masses used in the lattice remain above the physical values. The EM symmetry-breaking is 
of different nature as the mass, the former coming from vector-like interactions while the later being of scalar type. 
However, as we have discussed in the previous paragraphs, their symmetry breaking effects on certain observables are 
similar. Thus, the isovector part in (14) increases the masses of the charged mesons, according to ^ and ([7|, while 
the isovector and isoscalar both mix with the mass and contribute to {qq). There are other arguments pointing in 
the same direction when EM interactions are switched on. At finite temperature, the EM pure thermal corrections to 
the condensate also increase its absolute value for any temperature On the other hand, the condensate increases 
under the influence of an external magnetic field ei7, which can be also understood as the reduction of the free energy 
e ^ mq{qq) < (to leading order) needed to compensate for the EM energy increasing AesM ^ i^H)'^ /2 > [17]. 

Our purpose here will be to explore the consequences of that EM ferromagnetic behaviour to leading order. If the 
vacuum response is ferromagnetic, certain bounds for the EM LEG involved should be satisfied. We will derive those 
bounds and show that they are independent of the low-energy scale and thus can be checked in terms of physical 
quantities. Next we will check that the bounds are satisfied for the different estimates available for the EM LEG, with 
more detail for the SU{3) case in section IV D where we also discuss the gauge independence of our results. This will 
provide a consistency check for the ferromagnetic behaviour. 

An important comment is that we will discuss the ferromagnetic- like condition on the EM correction to (qq) and, 
as explained above, the splitting of the e = and e ^ parts in QGD-I-EM is ambiguous [39l[40]. This does not affect 
the low-energy representation of the condensates, which can be written in terms of physical quantities such as meson 
masses and decay constants. However, we will test our bounds with the EM LEG estimates obtained by matching 
low-energy results with the underlying theory |34fl37j . Therefore, we have to be consistent with the prescription 
for charge splitting followed in those works. This amounts to the direct separation of the e = part, which still 
may contain residual charge and fiQ QCD scale dependence through running parameters. The consequence is that 
the EM LEG thus defined are in general ^o-dependent, as discussed in [ID]- Therefore, those estimates are reliable 
only if there is a stability range where the dependence on /io is smooth and lies within the theoretical errors I34j . 
Actually, such stability range criteria are met for the LEG involved in our analysis (see section IV D I. Within that 
range, our identification of the dependent part in (qq) is consistent with the splitting scheme fol owed in those 
works. Actually, in that scheme Fq is /zo-independent and the e^-dependent part of the /xq running of Bomu,d is the 
same as that of m„,d in perturbative QGD-I-EM [30]. Hence, it is consistent to assume that the ChPT leading order 
of (qq) = —2BoF§ + . . . does not introduce any residual dependence when performing the charge splitting in the 
low-energy expression. Using a different splitting prescription would lead in general to different bounds and a different 
definition of the EM LEG. For instance, an alternative splitting procedure is introduced in |40| by matching running 
parameters of the theory with those of the = one at a given matching scale /ii. In that way, the EM part can 
be chosen as /io independent but it depends on the matching scale fii. We will not consider that splitting here, since 
there are no available theoretical estimates for the LEG defined with that scheme. The scale dependence for the LEG 
in either scheme is roughly expected to be logarithmic. 

Having the above considerations in mind, and going back to the case of physical quark charges, we separate the 
EM corrections to the condensate through the ratio: 
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= 1 + 2 - A..±] + e'JCUf^) + O (/) 



(15) 



with 



I'i defined in (|13k and where we have expanded in 5„ = ~M^o) /M^o — 0.1, as ^ij^^ ~ Mtt" — 

(^TT-M^o (j/^o/i^^) + 0{^), which is numerically reliable and can be performed in addition to the chiral expansion, 
in order to simplify the previous expression. 



We note that in 5/7(2) and to this order, the ratio (15 1 is not only finite and scale-independent but it is also inde- 



pendent of the not-EM LEG, including the contact /ii, /13, and therefore free of ambiguities related to the condensate 
definition. In fact, this ratio is also independent of Bq, unlike the individual quark condensates, which have only 
physical meaning and give rise to observables when multiplied by the appropriate quark masses, since m^iJo ^ Mf. 
In SU{2), the above ratio does not depend either on the mass difference md — niu, i.e., it depends only on the sum m 
and its deviations from unity are therefore purely of EM origin. All these properties make the ratio ( 15 ) a suitable 



quantity to isolate the EM effects on the condensate. Thus, the ferromagnetic-like nature of the chiral order parameter 
(qq), within its low-energy representation, would require that this ratio is greater than one, or equivalently to this 
order, d{qq) /de^ > 1. That condition leads to the following lower bound for the combination of EM LEG involved to 
this order, neglecting the C(5j) in (15) which changes very little the numerical results: 
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(16) 



We remark that the bound ( 16 1 is independent of the low-energy scale /i at which the LEG on the left-hand side 
are evaluated as long as the same scale is used on the right-hand side. Thus, it provides a well-defined low-energy 
prediction, expressed in terms of meson masses. The LEG on the left-hand side could be estimated by fitting low- 
energy processes or theoretically from the underlying theory, with all the related subtleties commented above. 

The condition (16) and the corresponding ones for SU{3) that will be derived in section IV D are obtained as a 



necessary condition that the LEG should satisfy if the QGD physical vacuum is ferromagnetic. This positivity condition 
on the quark condensate probes the vacuum by taking the mass derivatives ^ through the external source method so 
that the quark masses have to be kept different from zero and in that way the explicit symmetry breaking corrections 
are revealed in the condensate. If one is interested in the chiral limit, it must be taken only after differentiation, 
i.e., directly in eq.(15). In that case, it is not justified to perform the additional (5^ expansion in the charge because 
the e = masses vanish and one would be left only with the /i^± contribution in the r.h.s. of ( |15[ ), now with 
M^-t = 2Ce^/F^. That would actually give a larger negative value for the lower bound, coming from the smallness of 
the charged part of the pion mass (see details below) so that the bound in the chiral limit is less predictive. The fact 
that our bounds depend on quark masses is similar to other bounds on LEG obtained from QGD inequalities |4¥j . 

Nevertheless, the main physical interest is to test this bound for physical masses, using different estimates of the 
LEG in the literature. Thus, as a rough estimate, setting /i — Mp ~ 770 MeV and with the physical pion masses 
M^± ~ 139.57 MeV, M^o ~ 134.97 MeV, the bound ^ gives 5 [/cg(Mp) + fcg(Mp)] -Ffcf > -6.32(5.62) x 10"^ taking 
F = 87.1(92.4) MeV. This is a bit more restrictive than the "natural" lower bound -6.93 x 10"^ for the above LEG 
combinati on, o btained by setting all of them to — l/(167r^). The chiral limit gives —0.17 (with the value of C discussed 
87.1 MeV) i.e., much less restrictive, as commented above. More detailed numerical analysis 



m section 



IIA and F ■ 



will be done for SU (3) in section IV D 



On the other hand, the maximum value for the ratio ( 15 ) for the fc[ 



within "natural" values is obtained by setting 
the three of them to fc[(yLt = Mp) = l/(167r^), giving jl^j^^o = 1.0054, which gives an idea of the size of this correction. 
We remark that the term proportional to j/^ on the ratio ( 15 ) comes directly from the dependence of the pion masses 
on e^, so that it parametrizes the corrections in the condensate coming from any source of pion mass increasing, not 
only the EM one. Therefore, the same result can be used in order to provide a rough estimate of lattice errors in 
the condensate due to including heavier pion masses as lattice artifacts. In some lattice algorithms like the staggered 
fermion one, the situation is very similar to the mass differences induced by the charge terms. In that formalism, the 
finite lattice spacing induces terms |3S] that break explicitly the so called taste symmetry (four different quark species 
or "tastes" are introduced for every quark flavour) leaving a residual U{1) symmetry, pretty much in the same way as 
the charge term in ([2|. As a rough estimate, we can then replace 2Ce^/F^ by the corresponding 5-,^ from the lattice, 
obtained as the difference between the mass of the lightest lattice meson and the true pion mass. For a lattice pion 
mass of about 300 MeV, the v^^ term in (15) gives a correction of about 6%, which for a condensate value of (250 
MeV)'^ represents about (5 MeV)'^, which is within the order of magnitude quoted in 31]. Nevertheless, it should be 
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taken into account that the staggered ChPT ^3] has a much richer structure than the EM terms considered here and 
in particular there will be other operators contributing to the condensates at tree level, multiplied by the pertinent 
low-energy constants. If those constants are of natural size, we expect the size of the corrections to the condensate to 
remain within the range quoted above. 

IV. THREE FLAVOUR QUARK CONDENSATES. 
A. Results for light and strange condensates. 

In the SU{3) case, we derive to one loop in ChPT the light condensate {qq)i = {uu + dd) and the strange one (ss), 
taking into account both m„ — md and e 7^ corrections. Apart from the kaon and eta loops, an important distinctive 
feature in this case is the appearance of the ■n'^rj mixing term with the trec-levcl mixing angle e defined in (|6]) which is 
one of the sources of isospin-breaking corrections. The results we obtain for the condensates with all the corrections 
included are the following: 

(gg)f = {uu + ^ _2F2So |i + ^ [m {2Ll{^Ji) + H-^{^i)) + 4(2m + m,)Ll{pi)] + e^K.l+{^i) 

- ^ (3 - sin^e) ^^0 - 2/i^± - /iifo - - ^ (1 + sin^e) /i^ + O (17) 
{uu-ddf''^^^ - 2F■'B^{^{ma-m^){2Ll{^l)+Hl{^J))-e'lCl_{^l) 



sin 2e 



{ss) - -^25o^ + ^[m,(2L^(M) + i^2"(M)) + 4(2m + m,)Lg(A^)] + e2/C:(A^) 



[^ttO sin^ e + cos^ - 2 [^^± + ^iko] + O (p"^) | 



(18) 



(19) 



where we use the notation ( 13 ) and 



/CKm) = l[3{Kr + K^,{fi))+K^,(p)+Kl,{f,)]. (20) 

Note that in some of the above terms we have preferred, for simplicity, to leave the results in terms of quark instead 
of meson masses. An important difference between the SU{2) and SU{3) cases is that now there are loop corrections 
in {uu — dd) , where eta and pion loops enter through the mixing angle and kaon ones through the charged- neutral 
kaon mass difference. We have checked that the results are finite and scale-independent with the renormalization of 
the LEG given in Appendix [A| and that they agree with [3] for e = 0. Some unpublished results related to the SU{2) 
and SU{?)) isospin-breaking condensates can also be found in |50j . 



Numerical results for the condensates to this order can be found in [26] . As explained in section II A the effect of 



the Kl constants ( 20 1 in the condensates is not fully considered in that work, where the EM contributions are included 
through the corrections of Dashen's theorem |34j . so that only the Kl combinations entering mass renormalization 
appear. We will use then our results with all corrections included to estimate the range of sensitivity to the Kl of 
the condensates, analyzing the possible relevance for the fit in [26 . Our results are displayed in Table |l] As discussed 
in section III A we take the same input values Lg = 0, m^/m = 24 as in 26 as well as the assumption iJj — 2Lg, and 



the output values of B^m^ ^ si Wu/m^;, F, Lg from their main fit. In the second and third columns of Table |T] we give 
the results with all the EM Kl fixed to their minimum and maximum "natural" values. Since the Kl appear all with 
positive sign in ( 20 ) , the absolute values of the condensates obtained in this way are, respectively, lower and upper 
bounds within the natural range. We compare with the results quoted in |26j to the same C(p*) order (fourth column) 
for their main fit and we also show for comparison the results in the isospin limit e = 0, m„ = (fifth column). 
Our results agree reasonably with P6] . although we note that the values in that work lie outside the natural range 
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for the individual condensates. The largest relative corrections are about 2% for the light condensates and about 4% 
for the strange one. These isospin-breaking corrections are significant given the precision of the condensates quoted 
in [55]. On the other hand, the corrections lie within the error range quoted by lattice analysis [31]. In turn, note 
the bad ChPT convergence properties of the strange condensate, somehow expected since (ss) is much more sensitive 
to the strange quark mass ms than the light condensate [27] and therefore the large strange explicit chiral symmetry 
breaking nis is responsible in this case for the spoiling of the ChPT series, based on perturbative mass corrections. For 
the vacuum asymmetry — 1, the natural values band cover the result in [26| . although the numerical discrepancies 
in that case are relatively larger, between 15% and 24% for the lower and upper limit of the Kl respectively. Recall 
that this quantity vanishes to leading order in ChPT, according to (18 1, so that we expect it to be more sensitive to 
the Kl correction, which in this case comes mostly from the combination Kg + KIq. Nevertheless, it is worth noting 
that the results [26j imply (dd) / (uu) > 1 and {ss)/{uu) > 1, both in disagreement with many sum rule estimates of 
the condensate ratios [53] ■ Not surprisingly, we have the same discrepancy, since we use the same ChPT approach 
and the same numerical constants, except for the Kl corrections. The discrepancy in the relatively large value of 
{ss)/{uu) comes possibly from the bad convergence of the ChPT series for the strange condensate, which in addition 
is very sensitive to the choice of TJj EH]- The light condensates converge much better and although the sign of 
(dd) I (uu) — 1 is under debate, its magnitude is very small. In the latter case, our present calculation may become 
useful since the Kg + K^q contribution may change the sign of the vacuum asymmetry, although its precise value 
to fit a given prediction for {dd) / {uu) — 1 would still be subject to the H2 value. For this reason, it is important 
to make predictions for quantities which are independent of this ambiguity, as we have done in section |III| and as 
we will do in section |IVB[ where the sum rule for condensate ratios will allow to make a more reliable estimate of 
the vacuum asymmetry including both sources of isospin breaking. Finally, we comment on the numerical differences 
by considering the more recent low-energy fits in [28|. Still keeping H2 = 2Lg, these new values for the low-energy 
parameters increase considerably the total and strange condensates, which to ©(p**) give {qq) /{2BqF'^) ~ —2.15 and 
{ss) /{BoF"^) = —2.79. These higher values are mostly due to the much smaller F = 65 MeV, obtained in the main fit 
of [28] to accommodate a rather high L4 also with a large error L4 = (0.75 ± 0.75) x lO"'^ (an output result in ,28, ) . 
With the previous value F = 87.1 MeV but keeping the rest of LEC and masses as in [28] we get {qq) / {2BoF'^) ~ —1.63 
and {ss) / [BqF"^) = —1.99. The EM corrections remain of the same size and therefore their relative effect is somewhat 
smaller. As commented before, to„ 7^ isospin breaking is not implemented in those new fits and EM corrections 
are included only in kaon masses. 





-"-7-10 - - JE^ 


-"-7-10 - 


126] value Oip^) 


Isospin limit 


-{uu)o/{BoF'^) 


1.278 


1.292 


1.271 


1.290 


-{dd)o/{BQF'') 


1.297 


1.305 


1.284 


1.290 


-{ss}o/{BqF'') 


1.899 


1.907 


1.964 


1.904 


{dd) _ 
{uu) 


0.015 


0.010 


0.013 






TABLE I: Results for quark condensates. We compare with the values of (26j to 0{p'^) using the same set of low-energy 
parameters as in the main fit of that work, except the Kl , which we consider at their lower (second column) and upper (third 
column) "natural" values. We also quote the values in the isospin limit to the same chiral order. 



B. Sum rule corrections 



As noted in [3J , for to„ ^ one can combine the isospin-breaking condensates into a sum rule relating the isospin 
asymmetry {dd) / {uu) with the strange one {ss)/{uu). Such relation is phenomenologically interesting because it does 
not include contact terms and hence is suitable for numerical estimates on the size of the isospin-breaking corrections. 
Our purpose in this section is to discuss the EM e 7^ contribution to that sum rule. To leading order in rriu — rrid 
and we find: 



A 



SR ■ 



{dd) 



ma 



rrid 



1 



- TO 16tt^F^ 
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Ml 



Ml + Ml\og 



Ml 
Ml 
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{Kl{t,)+Kl,{y.)) 



(21) 
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The last term, proportional to is scale- independent and is the charge correction to the result in [3]. With the 
numerical set we have been using, the nid — rriu term in the right hand side gives —3.3 x 10~^, whereas the term 
gives -3.37 x lO^^ with K^{Mp) + Kl^iMp) = l/(87r2) and -9.4 x IQ-* with Kl{Mp) + Kl^iMp) = 2.7 x 10-^ the 
central value given in Therefore, the charge term above is of the same order as the pure QCD isospin correction 
and must be included when estimating the relative size of condensates through this sum rule. In fact, using the values 
quoted in ^29j mu/md ~ 0.55, ms/md — 18.9 and {ss)/{uu) — 0.66, we get from (21) with physical pion and kaon 
masses: 



-0.015 < 



(dd) 



1 < -0.009, 



where the lower (upper) bound corresponds to the natural value Kg + KIq — -|-(— )l/(87r^), while the value without 
considering the charge correction is —0.012 and the value quoted in [53] collecting various estimates in the literature is 
—0.009. The inclusion of the charge corrections may then help to reconcile this sum rule with the different condensate 
estimates available. In fact, through this sum rule we see that ChPT is also compatible with the asymmetries (dd) / (uu) 
and (ss) / (uu) being both smaller than one (see our comments in section IV A). Note that the ferromagnetic- like 
arguments used in sections III and IV D cannot be applied to {uu — dd) , which does not behave as an order parameter 
under chiral transformations, since it is not invariant under SUv{'2). Finally, we recall that estimates based on the 



sum rule (21) are more precise than the ones we have made directly from the condensates in section IV A since this 



sum rule is free of the H2 ambiguity. 



C. Matching of low-energy constants 



Our aim in this section is to explore the consequences of including the two sources of isospin breaking for the 
matching of the LEG involved in the condensates. For that purpose, we perform a l/rris expansion in the SU{S) sum 
and difference condensates given in (17)-(18). Matching the 0(1) and 0{\ogms) terms with the corresponding SU{2) 
expressions in (lO)-(ll) yields the following relations between the LEG, for the sum and difference of condensates 
respectively: 



2ACo [ll{^i) + hlifi)]+e'F'mt,) 



2Mi 



16L^(M)+4Lg(M) + 2i72^(M) 



1^3 



18 



2 



(22) 



Bo{md - TO„)/l3 — fcy 



Bo{md - TO„) 
2e^F^ ' 



4Llif,) + 2H^,ifi) 
3C 



2 



2F^ 



Vko 



967r2 



(23) 



In the above expressions, we have displayed the SU{2) contribution on the left hand side and the SU{Z) ones 
on the right hand side, with IC^il^) and /CajJ/i) given in (12) and (20). Note that the 1/ms expansion has been 
implemented also in the tree level relations ([7|, so that M^o = (w„ + md)BQ + 0{l/ms), M^o = B^rag -\- 0(1) and 
A/^ — ABQmg/'i + 0{l). It is important to pomt out that the pion mass charge difference is not negligible in the l/m^ 
expansion, and for that reason we keep M^o in (22 1. For kaons, it is justified to consider the charge contribution 



negligible against the dominant term, so that at this order M^± and M^o are not distinguishable. 

In the sum matching relation (22), the isospin corrections are not very significant. The mass difference m„ — md 



does not appear in the neutral and kaon masses to leading order in l/mg and the charge correction, although being 
of the same chiral order as the Af^o term, numerically e^F^/M^o ~ Ce^/(F^Af^o) — 0.05. However, in the difference 
matching (23 ), the to„— corrections contribute on the same footing as the EM ones and are numerically comparable. 



The above matching relations can be used directly for the approximated LEG (estimated theoretically or fitted to 
data) and for physical masses, since the difference with the tree level masses and LEG is hidden in higher orders. 
On the other hand, for the tree-level LEG, i.e. the GhPT 0{p^) lagrangian parameters, since they are formally 
i nde p end ent of the light quark masses, we can just take the chiral limit m„ — md = in the abo ve expre ssions 
(22)-(23l and read off the corresponding matching of the contributions. Using the latter again in (22)-(23) gives 
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then independent relations between the tree-level LEG involved at = and the EM ones. Doing so, the EM and 
not-EM LEG combinations decouple and the results are compatible with those obtained in [3. for e = and in 03] 
for e ^ (setting to„ — = from the very beginning): 



5(fc5^(/i) + fcg(M)) 



16L^(m) + 4i^(M) + 2H;{^,) - ^ - 1^ 



2 



1 



967r2 



3C 



6{K, + K]^{p)) + A{K^,i^) + - j^i^KO 

3C 



(24) 



where the v functions are evaluated exactly in the chiral limit, i.e, for Mj^o — BquIs 
and third equation coming from ( 22 ) and the second and fourth from ( 23 1 . 



and = 4B0TO5/3, the first 



Our first conclusion is then that to this order of approximation, the formal matching of the condensates is consistent 
with the matching relations previously obtained. In other words, mass and charge terms can be separately matched. 
This would be no longer true at higher orders where for instance e^irriu — rnj) contributions may appear. 

Although the relations (22) and (231 reduce to (24) in the chiral limit for the tree-level LEG, it is better justified to 



use the original expressions (22)-(23l when dealing with physical meson masses and when the LEG are obtained either 
from phenomenological or theoretical analysis. The LEG obtained in that way are approximations to the lagrangian 
values and consequently they depend on mass scales characteristic of the approximation method used. For instance, 
the LEG obtained by phenomenological fits are sensitive to variations both in rh and in m„ — vrid |26| . in resonance 
saturation approaches they depend on vector meson masses [371 138j which themselves depend on quark masses and 
in the NJL model some LEG like K^^ depend on the scale where the quark masses are renormalized |31|. We do not 



expect large differences between using the general matching relation ( 22 ) or the first and third equations in ( 24 1 , since 
the latter can be understood also as the e = limit of the former and we have seen that this is numerically a good 
approximation. However, that is not so clear for (22 1 where the two isospin-breaking contributions are of the same 
order, both in the chiral expansion and numerically. 



Finally, we can use the previous matching relations to estimate numerically the SU{2) condensates in (lOl-(ll) 
without having to appeal to the values of the SU{2) LEG. Doing so we obtain -{uu + dd)^^^'^^ / BqF'^ ~ (2.16, 2AE) 
and (uu — dd)^"'''^^/BoF^ ~ (0.014,0.02) where we indicate in brackets the natural range of the EM LEG, to be 
compared to -{uu + dd)^^^^y B^^F"^ ~ (2.58,2.6) and {uu - dd)^^^ B^^F^ ~ (0.013,0.018) from Table|l] The larger 
difference in {uu + dd) comes from the 0{ms) and 0{ms log rris) terms in the I/to^ expansion, which were separated 
when doing the matching and which are absent in the condensate difference. In fact, the numerical contribution of 
those terms to — {uu + dd)^^''^^ / BqF"^ is about 0.41, which explains perfectly the numerical differences and confirms 
the idea that in standard GhPT the light condensates calculated either in the SU{2) or in the SU (3) cases give almost 
the same answer near the chiral limit. This may be not the case in other scenarios of chiral symmetry breaking |32j . 



D. EM corrections and 5(7(3) LEG bounds 



We have seen in the SU{2) case that the EM ratio given in ( 15 ) is a relevant physical quantity allowing to establish 
a constraint for the EM LEG based on explicit chiral symmetry breaking. The same argument applied to the 5'C/(3) 
case leads also to a constraint on the EM LEG obtained from the full condensate {qq) = {uu + dd+ss), which behaves 
as an order parameter, being an isosinglet under 5'J7y(3). In addition, we can still consider the light condensate {qq)i 
as the order parameter of chiral transformations of the SU (2) subgroup, which in principle will lead to a different 



constraint. In fact, the latter is nothing but the constraint obtained in the SU{2) case (16 1, once the equivalence 
between the LEG obtained in section |IV G| is used. As for the full condensate, it should be kept in mind that the 
large violations of chiral symmetry due to the strange quark mass may spoil our simple description of small explicit 
breaking. As commented above, this reflects in the large NLO contributions to the strange condensate, which in 
the standard GhPT framework depends strongly on m^, unlike the light condensate. Therefore, the bounds of the 
LEG obtained for the full condensate are less trustable, since neglecting higher orders, say of ©(e^m^) is not so well 



justified for (ss) . Proceeding then as in section III where the same prescription of charge splitting when comparing 



with QGD approaches is understood, we calculate the ratios: 
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S(7(3) 



SU{3) 



l + ^[6{Kj + Kliii)) + 5 {Klin) + KM)] 



1 + -e^ [2{Kl{^i) + K{^{^,)) + 3(/w + if8(M))] 



(25) 



(26) 



0.008, which 



where, as in the SU{2) case, we have expanded in Sjr and also in Sk 

allows to express the results in terms of the full tt* and masses. Otherwise we should take into account that now 
M^±{e = 0) 7^ Mt^o, unlike the SU{2) case, and Afjc±(e = 0) 7^ Mj^o, by terms of order — m„. This is important 
when using this result for numerical estimates, since, as discussed before, the separation of the e — correction to 



the masses is formally not unique. As in SU{2), the ratios (25)-(26) are finite and independent of the scale /i, of Bq 
and of the not-EM LEG, so they are free of contact ambiguities. 

As in section we want to explore the consequences of the ferromagnetic nature of the physical QCD vacuum 
under explicit chiral symmetry breaking for the EM LEG. Here also the charge-mass crossed terms in the fourth 
order lagrangian (A2| give explicit breaking contributions to the quark condensate coming now from the isoscalar, 
isovector and strangeness part of the charge matrix. For physical quark charges, demanding that the ratios (25l-(26) 



are greater than one we obtain the following EM bounds, to leading order in the chiral expansion and in Sj^tSk- 
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{qq)i ^ 6 {Kj + KliiJi)) + 5 {Kl{^,) + if[o(M)) > — (2j.,± + 

•in 

(qq) 2{Kj + Kli^i)) + 3 (i^9^(M) + if[o(M)) > (^.± + ^K±) 



(27) 
(28) 



We remark that these constraints are independent of the low-energy scale /j,. It is also clear that the light bound 



(27) is nothing but the one obtained in the SU(2) case (16) once the equivalence between the LEG given in the third 
equation of (24) is used. In Figure [T] we have plotted these two constraints in the (ify -I- K^) — {Kg + KIq) plane at 
fi = Alp and within the natural region. We have used the same numerical values for the tree level LO masses, C and 
F as in previous sections. We observe that the bound on the full condensate is more restrictive than the light one in 
that range. However, as we have commented above, it is also less trustable, due to the large distortion of the chiral 
invariant vacuum due to the strange mass. Both bounds give also a more restrictive condition than just the natural 
size. 
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FIG. 1: Regions in the LEG space constrained by the bounds on the light condensate (271, above the full blue line, and the 



full one (28), above the dashed red line. The LEG are renormalized at ^ = Mp and are plotted within the natural range. 



Let US now check these bounds against some estimates of the K^ in the literature. We start with the sum rule 
approach for K^ in j37]. In that work, K^ = K^{Mp) = 0, but what is more relevant for us is that the combination 
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KIq at any scale is gauge-independent despite being both Kg, K{q 



dependent on the gauge parameter 5, as one 
can readily check from the explicit expressions given in I37[ (c.f. their eqns. (94) and (95)). This is an interesting 
consistency check of our present bounds (27 1 and (28), which are gauge independent in addition to their /i low- 
energy scale independency commented previously, supporting their validity and predictive character. Numerically, 
the constant could not be estimated in [37| due to the slow convergence of the integrals involved, but they provide 
a numerical estimate for if[g(Mp) — 5.2 x 10~^ at ~ 0.7 GcV and ^ = 1, for which we get Kl{Mp) > —0.021 from 
(271 and Kg{Mp) > —0.015 from (28 1. See our comments about the /io scale dependence in section 

In [35 resonance saturation gives iQ^Mp) = 0, which is compatible with our present bound, 
incompatible with a previous value for K^{Mp) = (— 4±1.7) x 10~^ obtained in [3]. The possible reasons to explain this 



III 

This is apparently 



difference were discussed in 



That value for Kg is compatible with our bounds as long as 6KT+5{Kg+KlQ) > —0.05 
0.02 from pSl). 



from ^ and + 3(iCJ _ 

In 1^ based on large Nc and the NJL-model, the LEG estimates give = 0, Ks{Mp) = (-0.8 ± 2.0) x 10"^ {Kr 
and Ks are 0{1/Nc) suppressed) and K'g{Mp) + Kl^{Mp) = (2.7 ± 1.0) x lO'^, aU of them at = 0.7 GeV. These 
values are also compatible with both bounds ( 27 ) and ( 28 ) . We note that in |36] , where the short-distance contributions 
are evaluated as in 



the explicit expressions given for the LEG show again that K-i,K^ and + K^q are gauge 



independent. Furthermore, 
in the /iq scale around /io = 
bounds, the comparison with those works is robust concerning the gauge and QGD scale dependence. 



Kq and KIq, dominant for large Nc in that approach, show a rather large stability range 
: 0.7 GeV [34j[36]. Since Kg + KIq is the only combination surviving for large Nc in our 



V. CONCLUSIONS 



In this work we have carried out an analysis of strong and electromagnetic isospin-breaking corrections to the 
quark condensates in standard one-loop GhPT, providing their explicit expressions and studying some of their main 
phenomenological consequences for two and three light flavours. 

Our results have allowed us to analyze the sensitivity of recent isospin-breaking numerical analysis of the condensates 
to considering all the EM LEG involved. The effect of those LEG is smaller for individual condensates than for the 
vacuum asymmetry, where they show up already in the leading order. These corrections lie within the error range 
quoted in lattice analysis. Our analysis can also be used to estimate corrections to the quark condensate coming from 
lattice artiflcially large meson masses. 

We have shown that if EM explicit chiral symmetry breaking induces a ferromagnetic-like response of the physical 
QGD vacuum, as in the case of quark masses, one obtains useful constraints as lower bounds for certain combinations 
of the EM LEG, both in the two and three flavour sectors. We have explored the consequences of this behaviour for the 
ratios of e 7^ to e = light and total quark condensates, which are free of contact-term ambiguities, and for a given 
convention of charge separation. The large GhPT corrections to the strange condensate make the constraints on the 
full condensate less reliable. In this context, we have discussed the different sources for EM explicit chiral symmetry 
breaking and isospin breaking terms, by considering formally arbitrary quark charges. Thus, there are chiral symmetry 
breaking terms proportional to the sum of charges squared, coming from crossed charge-mass contributions in the 
effective action, which show up in the vacuum expectation value. In accordance with the external source method, 
we keep the quark masses different from zero to account correctly for all the explicit symmetry breaking sources. 
The chiral limit can be taken at the end of the calculation. The bounds obtained are explicitly independent of the 
low-energy scale /i, providing then a complete and model-independent prediction at low energies. However, when this 
low-energy representation is compared with theoretical estimates based on QGD, one has to take into account that 
due to the convention used in the charge separation, the estimated LEG depend on the QGD renormalization scale 
/io, as well as being gauge dependent. Our bounds are numerically compatible with those estimates, based on sum 
rules, resonance saturation and QGD-like models, within the stability range of /xq where those approaches are reliable. 
Furthermore, the LEG combinations appearing in our bounds are gauge independent. We believe that our results can 
be useful in view of the few estimates of the EM LEG in the literature. 

We have found that the EM correction to the sum rule relating condensate ratios is of the same order as the 
previously calculated e = one, and therefore must be taken into account when using this sum rule to estimate the 
relative size of quark condensates. We have actually showed that using the complete sum rule, which is also free of 
contact terms, yields a GhPT model-independent prediction for the vacuum asymmetry compatible with the results 
quoted in the literature. 

Finally, we have performed a matching between the SU(2) and SU{S) condensates, including all isospin-breaking 
terms. Matching the sum and difference of light condensates gives rise to matching relations between the LEG 
involved, where EM and not-EM LEG enter on the same footing in the chiral expansion. These matching relations 
may be useful when working with physical masses and LEG estimated by different approximation methods. In the 
case of the sum, the charge contribution is numerically small with respect to the pion mass one, but in the difference 



15 



the two sources of isospin-breaking are comparable. Taking the chiral hmit, EM and not-EM constants decouple and 
the matching conditions are compatible with previous works for the LEG in the lagrangian, which are defined in this 
limit. 

Appendix A: Fourth order lagrangians and renormalization of the LEG 

We collect here some results available in the literature and needed in the main text. To calculate the quark 
condensates to NLO one needs the £p4_|_p2g2_|.g4 lagrangians to absorb the divergences coming from loops with vertices 
from £p2_|_g2. We denote by a superscript qq the relevant terms in the lagrangian, which are those containing the 
quark mass matrix. For SU{2) they are [5, 6]: 

^l' = + + Je^i + '^3)tr[x'] + \{h^- h,) det(x), 



C%2 = (kMxiU + U^)MQ^] + ketr[x{U + U^)MQUQU^] + kM{xU^ + Ux)Q + {xU + Uh)QMQ] 

where x = '^BqM, whereas for SU(3) [4]: 

£«I =LetY[x{U + C/t)]2 + LstvixUxU + xUhU^] + H^x^] 



[KjiT[x{U + U^)]tr[Q^] + KsiT[x{U + U^)]ir[QUQU^] + KQir[{xU + U^x + + Ux) (A2) 



i^iotr [{xU + U^x) QU^QU + {xU^ + Ux) QUQU^] j , 
and — for both cases. 

In order to renormalize the meson loops it is necessary to separate the low energy constants appearing in the NNLO 
lagrangian in finite and divergent parts. The renormalization of the LEG involved in the calculation of the SU(2) 
condensates is given by [21 13 E] : 

hi =hl{fi) + 5iX, 
ki + CTjA, 



with 73 = — i, (5i = 2, = 0, and 0-5 = — | — |Z, cg = j + 2Z and 0-7 = 0, for physical quark charges e„ — 2e/3, 
ed = -e/3, being Z := ;pr. 
expressed as 



= —e/3, being Z := -S^. The part that diverges in d = 4 dimensions is isolated from the loop integrals and is 



where — r'(l) is the Euler constant. As for the SU(3) ones, we have [31 13] 

U =LI(M)+r.A, 

=iJ[(M) + A,A, 

with Te = Fg = ^, A2 = ^, and Ey = 0, Eg = Eg = -\, Eio = | + fZ. 
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